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Abstract 

We consider a strong Markov process with killing and prove an approximation 
method for the distribution of the process conditioned not to be killed when it is 
observed. The method is based on a Fleming- Viot type particle system with rebirths, 
whose particles evolve as independent copies of the original strong Markov process 
and jump onto each others instead of being killed. Our only assumption is that 
the number of rebirths of the Fleming- Viot type system doesn't explode in finite 
time almost surely and that the survival probability of the original process remains 
positive in finite time. The approximation method generalizes previous results and 
comes with a speed of convergence. 

A criterion for the non-explosion of the number of rebirths is also provided for 
general systems of time and environment dependent diffusion particles. This in- 
cludes, but is not limited to, the case of the Fleming- Viot type system of the ap- 
proximation method. The proof of the non-explosion criterion uses an original non- 
attainability of (0,0) result for pair of non-negative semi-martingales with positive 
jumps. 

1 Introduction 

Markov processes with killing are Markov processes which stop to evolve after a random 
time, called the killing time. While the behavior of a Markov process after this time is 
trivial, its distribution before its killing time represents a substantial information and 
thus several studies concentrate on the distribution of a process conditioned not to be 
killed when it is observed (we refer the reader to the extensive bibliography updated by 
Pollett (2o| . where several studies on models with killing mechanisms are listed). The 
main motivation of the present paper is to provide a general approximation method for 
this distribution. 

There are mainly two ways of killing a Markov process, both of them being handled 
in the present paper. The first kind of killing occurs when the process reaches a given 
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set. For instance, a demographic model is stopped when the population becomes extinct, 
that is when the size of the population hits 0. The second kind of killing occurs after 
an exponential time. For example, a chemical particle usually disappears by reacting 
with another one after an exponential time, whose rate depends on the concentration of 
reactant in the medium. If the killing time is given by the time at which the process 
reaches a set, we call it a hard killing time. If it is given by an exponential clock, we call 
it a soft killing time. In order to formally describe a process with killing in a general 
setting, we consider a continuous time strong Markov process Z evolving in a measurable 
topological state space F U {d}, where d ^ F is absorbing. This means that Z evolves 
in F until it reaches d and then remains in d forever. The killing time of Z is defined as 
the hitting time of d and is denoted by tq = inf{t > 0,Zt = d}. This general setting can 
be used to model both hard and soft killing types. 

The main difficulty in approximating the distribution of a Markov process before its 
killing is that the probability of the event "the process is still not killed at time t " vanishes 
when t goes to infinity. Indeed, this is a rare event simulation problem and, in particular, 
naive Monte-Carlo methods are not well-suited (the number of still not killed trajectories 
vanishing with time, a Monte-Carlo algorithm would undergo an unavoidable discrepancy 
in the long term). Our approximation method is based on a quite natural modification 
of the Monte-Carlo method, where killed particles are reintroduced through a rebirth in 
F (the state space of the still not killed particles) in order to keep a constant number of 
meaningful particles through time. This modification, which appears to be a Fleming- 
Viot type particle system, has been introduced by Burdzy, Holyst, Ingermann and March 
in Q and is described in the following definition. 

Definition. Let Z he a strong Markov process with killing evolving in F U {d}. Given 
N > 2 and (xi, . . . ,xn) G F^ , we define the Fleming- Viot type particle system with 
N particles and evolving as Z between their rebirths, starting from {xi, . . . ,xn) and 
denoted by FV^(xi, . . . ,xj\[), as follows. The particles of the system start from 
(xi, ... ,3; at) at time and evolve as independent copies of Z until one of them 
is killed. Then 

- if two or more particles are simultaneously killed at this moment, or if a 
particle jumps simultaneously without being killed (so that a particle is killed 
and another one jumps), then we stop the definition of the particle system 
and say that the process undergoes a failure. 

- otherwise {i.e. if one and only one particle is killed and the other particles do 
not jump at this killing time) then the unique killed particle is taken from the 
absorbing point d and is instantaneously placed at the position of a particle 
chosen uniformly between the — 1 remaining ones; in this situation we say 
that the particle undergoes a rebirth. 

After this operation and in the case of a rebirth, each of the A^ particles lies in F. 
They then evolve as A^ independent copies of Z, until one of them is killed and so 
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on. This procedure defines the particle system FV^ {xi, . . . ,xj\f) in an incremental 
way, rebirth after rebirth and until it undergoes a failure. 



One particularity of this method is that the number of rebirths faces a risk of explosion 
in finite time: in some cases, that is for some choices of Z and N > 2, the number of 
rebirths in the particle system explodes to +00 in finite time. In such a case, the particle 
system is well defined up to the explosion time, but there is no natural way to extend 
its definition after this time (we refer the reader to i for a non-trivial case where the 
number of rebirths explodes in finite time with positive probability). This leads us to 
the following assumption (see also Section |3l where we provide a sufficient condition for 
a diffusion process to satisfy Hypothesis A(A^)). 

Hypothesis A(A^), > 2. A Markov process with killing Z is said to fulfil Hypoth- 
esis A(A^) if and only if the Fleming-Viot type particle system with N particles 
evolving as Z between their rebirths undergoes no failure and if its number of re- 
births remains finite in finite time almost surely, for any initial distribution of the 
particle system. 



In Section [21 we fix T > and consider a continuous time Markov process with killing Z. 
We state and prove an approximation method for the distribution of Zt conditionally 
to T < TQ. This method, based on Fleming-Viot type particle systems, does not require 
that Z satisfies Hypothesis A(N) for some N > 2. Instead, we assume that there exists 
a sequence {Z^)]^->2 of Markov processes such that Z^ converges in law to Zt when 
A^ — )• 00 and such that Z^ satisfies Hypothesis A(A^), for all N > 2. In particular the 
Fleming-Viot type system with N particles built over Z^ is well defined for any N > 2 
and any initial distribution. This particle system, denoted by FV^, evolves in and 
we denote by fi^ its empirical distribution at time t > 0. Our main result states that, if 
the sequence of initial empirical distributions {fJ.Q)N>2 converges to a probability measure 
fiQ on F, then the sequence of random probability measures {ij^t)^>0 converges in law 
to the distribution of Zt, initially distributed as and conditioned to T < rg. Also, 
we provide a speed of convergence for the approximation method. This result generalises 



the approximation method proved by Grigorescu and Kang in for standard multi- 



dimensional Brownian motion, by Villemonais |24| for Brownian motions with drift and 
by Del Moral and Miclo for softly killed Markov processes (see 0| and references therein) 
and by Ferrari and Marie j^, which regards continuous time Markov chains in discrete 
spaces. 

Our result gives a new insight on this approximation method through three key points of 
generalisation. Firstly, this new result allows both hard and soft killings, which is a natu- 
ral setting in applications: typically, species can disappear because of a lack of newborns 
(which corresponds to a hard killing at 0) or because of a brutal natural catastrophe 
(which typically happens following an exponential time). Secondly, we implicitly allow 
time and environment dependencies, which is quite natural in applications where indi- 
vidual paths can be influenced by external stochastic factors (as the changing weather). 
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Last but not least, for any > 2, we do not require that Z satisfies Hypothesis A(A^). 
Indeed, we only require that there exists an approximating sequence (-Z^)Ar>2 such that 
fulfils Hypothesis A(A^), VA^ > 2. This is of first importance, since the non-explosion 
of the number of rebirths is a quite hard problem which remains open in several sit- 
uations. It is for instance the case for diffusion processes with unbounded drift, for 
diffusions killed at the boundary of a non-regular domain and for Markov processes with 
unbounded rates of killing. Thus in our case, the three irregularities can be handled 
by successive approximations of the coefficients, domain and rate of killing. For these 
reasons, the approximation result proved in the present paper can be used as a general 
approximation method for the distribution of conditioned Markov processes. 



The approximation method being proved in this very general setting, the only remaining 
difficulty is to provide a sufficient criterion for Hypothesis A{N) to hold, for any N >2. 
In Section 121 we consider a multi-dimensional diffusion process Z and provide a sufficient 
criterion for Hypothesis A(A^) to be fulfilled for FV^ ^ the Fleming- Viot type system with 
A'^ particles evolving as Z between their rebirths. We allow both soft and hard killings 
and we allow the coefficients of Z to be time-inhomogeneous and environment dependent. 
This criterion generalises similar non-explosion and non-failure results recently proved 
by Lobus in 18 1 and by Bienek, Burdzy and Finch in j^l for Brownian particles killed at 



the boundary of an open set of R , by Grigorescu and Kang in |l3| for time-homogeneous 



particles driven by a stochastic equation with regular coefficients killed at the boundary 



of an open set and by Villemonais in |24l | for Brownian particles with drift. Other models 
of diffusions with rebirths from a boundary have also been introduced in iL with a 
continuity condition on the jump measure that isn't fulffiled in our case, in Il2l. where 



fine properties of a Brownian motion with rebirth have been established, and in [15[, [16|, 
where Kolb and Wiikber have studied the spectral properties of similar models. 

Our non-explosion result is a generalization of the previously cited ones and is actually not 
restricted to Fleming- Viot type particle systems. The ffist generalization axis concerns 
the Markov processes that drive the particles between their rebirths: our criterion allows 
a different state space and a different dynamic for each particle. Moreover these dynamics 
can be time-inhomogeneous environment-dependent diffiision processes, while previous 
results are restricted to time-homogeneous diffusion processes. The second aspect of the 
generalization is related to the rebirths mechanisms: we allow the rebirth position of a 
particle to be chosen with a large degree of freedom (in particular, not necessarily on 
the position of another particle) and the whole particle system is allowed to undergo 
a rebirth when a particle is killed. This setting includes (but is not limited to) the 
Fleming- Viot type particle system case, so that our non-explosion result validates the 
approximation method described above for time/environment dependent diffusions with 
hard and soft killing. For non Fleming- Viot type rebirths mechanisms, it remains an 
open problem to determine the existence and, if appropriate, the value of the limiting 
empirical distribution. 

The proof of the non-explosion is partly based on an original non-attainability of (0,0) 
result for semi-martingales, which is stated in Section |4] of this paper. Note that this 
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result answers a different problematic from the rest of the paper and has its own inter- 
est. Indeed, inspired by Delarue our non-attainability criterion generalizes existing 
non-attainability criteria by considering general semi-martingales (i.e. not necessarily 
obtained as solutions of stochastic differential equations). This essential improvement is 
a sine qua non condition for the development of the non-explosion criterion of Section 3. 

2 Approximation of a Markov process conditioned not to 
be killed 

We consider a cadlag strong Markov process Z evolving in a state space F U {d}, where 
9 ^ -F is an absorbing point for Z, also called the cemetary point for Z (this means that 
the process remains in {d} forever after reaching it). When Z reaches d, we say that it is 
killed and we set Tg = inf{t > 0, Zf = d}. In this section, we assume that we are given a 
random probability measure /io on F and a deterministic time T > 0. We are interested 
in the approximation of P^g {Zt £ ■ \ T < tq), i.e. of the distribution of Zt with initial 
distribution fiQ and conditioned to T < tq. 

Let {^j,q)n>2 be a sequence of empirical distributions with points which converges 
weakly to /ig (by an empirical distribution with N points fiQ , we mean that there exists 
a random vector (xi, ... ,2; at) in F^ such that fiQ has the same law as X^^Li ^Xi )- We 
assume that there exists a sequence {Z^)n>2 of strong Markov processes with killing 
evolving in FU {d}, such that Z^ fulfils Hypothesis A(A^) for any N >2 and such that, 
for any bounded continuous function / : F — )• M, 

E^. [f{Z^)l^^^.) K,, {f{ZT)lT<re) , (1) 

where Tq = inf{t > 0, Z^ = d} denotes the killing time of Z'^ . 

Remark 1. A particular case where ([1]) holds true (we still assume that fj,o is a random 
measure) is when is the empirical distribution of a vector distributed with respect to 
/i® ^ for any N > 2 and when, for all x G and all continuous and bounded function 
/:F^IR+, 

Prfi^) T— ^ Prfix). (2) 

Indeed, we have 

1=1 

where {xi)i>i is an iid sequence of random variables with law fiQ. By the law of large 
numbers, the first right term converges to almost surely and, by ([2]) and by dominated 
convergence, the second right term converges almost surely to no {Pt/)- 
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We define now the A^-particles system = as the Fleming-Viot type 

system whose N particles evolve as independent copies of between their rebirths and 



FV^(xi, ...,xn), where (xi, 



-.XN) 



with initial empirical distribution /^q (thus X 

is such that ^uq = ^/^^^=i^Xi)- Since we assumed that fulfils Hypothesis A(A^), 
the particle system is well defined at any time t > 0. Also, for any i G {1, . . . ,A^}, 
the particle X^'^ is a cadlag process evolving in F. We denote by iif the empirical 
distribution of Xi^ at time t > 0: 



1 ^ 

— > 5^i,f 
N ^ 



G Aii{F) almost surely. 



i=l 



where Aii{F) denotes the set of probability measures on F. 

Theorem 1. Assume that the survival probability of Z at time T starting with distri- 
bution fiQ is strictly positive almost surely, which means that P^,j (T < tq) > almost 
surely. Then, for any continuous and bounded function f : F ^ 11+, 



^i'^{f)-^^,Af{ZT)\T<T3). 



Moreover, for any bounded measurable function f : F ^ 11+, we have 



E 



l^^if) 



E^.{f{Z^)\T<Ti 



< 



2(1 + ^/2)11/11, 



N 



Remark 2. In [23j, Rousset consider a process Z which is only subject to soft killings, with 
a uniformly bounded killing rate. In this particular case, the author proves a uniform 
rate of convergence over all times T, using the stability of the underlying Feynman-Kac 
semi-group. Our result shall be used in further work to extend Rousset's results to many 
process with sufficiently stable associated Feynman-Kac semi-group. 

Remark 3. We emphasize that our approximation method only concerns the conditional 
distribution of the process. For results on the pathwise behaviour of the process condi- 
tioned not to be killed in the future, we refer the reader to the theory of Q-processes, 
which are processes evolving as Z conditioned to never be killed (see for instance |l7|). 

Remark 4- The quasi-limiting distribution, i.e. limf_>_|_oo P^jq {^t & ■\T < tq), is the 
subject of several studies and is closely related (in fact equivalent) to the concept of quasi- 
stationary distribution (see for instance (l9|). Theorem [1] implies that the quasi-limiting 
distribution, if it exists, is equal to limt-^^oo^^^N-^oo ^J't^ ■ Whether this limit is also 
equal to limN^ao^^T^t^+oo fJ-t^ (which has practical interests in numerical computation) 
remains an open problem in many situations and has been resolved quite recently in some 
particular cases (see for instance (5|,B,|i3 and 0)- 
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Proof of Theorem[J\ For any N > 2 and any i G {!,••• ,N}, we denote the increasing 
sequence of rebirth times of the particle X*'^ by 

For any t > 0, we denote by Al'^ the number of rebirths of the particle X'''^ before time 
t: 

4^ = max{n G N, < t}. 

We also set = X^^^ ^I'^ the total number of rebirths in the A^-particles system 
before time t. Since is assumed to fulfil Hypothesis A(A^), we have lim^^oo Tn^ = 
+00 almost surely and < +oo for all t > almost surely. 

For any N > 2, we denote by (P/^)t>o the sub-Markov semi-group associated with the 
Markov process Z'^ . It is defined, for any bounded measurable function / on F U {d} 
and any t > 0, by 

Pff(x) = (/(Zf , Vx E F U {d}. 

The proof of Theorem [1] is divided into three steps. In a first step, we fix > 2 and 
we prove that, for any bounded and measurable function / on F, there exists a local 
martingale such that 



N ^ 

N J 



=1 n=l j^i 



In a second step, we define the non-negative measure f/^ on F by 



Vi G [0,T]. 
(3) 



The measure is obtained by multiplying the empirical measure fi^ by (A^ — 1)/A^ at 
each rebirth. This loss of mass is introduced in order to compensate the last right term 
in ([3|): we prove that f/^ {f)~ fJ'o {^t^ f) is a local martingale whose norm is bounded 
by {1 + V2)\\f\\oo/VN, thus 



.N/j^\ _ ^1 ..N I 



E - l^^iP^ff) < (4) 

In step 3, we note that and //^ are almost surely proportional and conclude the proof 
of Theorem [1] by renormalizing and {^t ") iiiSl 

Step 1 : proof of decomposition ([3)) . 

Fix N > 2 and let / : F U {5} i— )■ B,-|_ be a measurable bounded function such that 
f{d) = 0. We define, for all t G [0,r] and z G F, 
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N 



The process is a martingale with respect to the natural filtration of Z 

and is equal to at time Tq on the event {tq < T}. Indeed, for all s,t > such that 
s + t <T, the Markov property implies 

and Pt^f{d) = for any t > yields 



Fix i G {!,..., iV}. For any n > 0, we define the process ^Mj'""^ - which implicitely 

depends on iV - by 

Mr = h^r^-MXl'"") - ^,,^.-(^:;^,.), with T^'^ := 0. 

Since X^'^ evolves as in the time interval [Tn^,T^'^i[, Mj'" is a martingale with 
respect to the natural filtration of the whole particle system and 

{0, if t < I.e. Ai'^ < n, 

MXl''') - if t G ^.e. Ai'^ = n, 

-'^p ^,N{X\^J^),l^t>T^^-^ t.e.n<A^' . 

Summing over all rebirth times up to time t, we get 

= Mx^'"") + E Mr + E v,.-(x;f^). (5) 

n=0 n=l 

Defining the local martingales 

]Mj=5^]Mrand Mt = -^M] 

n=0 1=1 

and summing over i G {!,..., A^}, we obtain 



i=l n=l 

At each rebirth time Tn^, the rebirth position of the particle X*'^ is chosen with respect 
to the empirical measure of the other particles. The expectation of -0 i,N{X^'^j^) condi- 
tionally to the position of the other particles at the rebirth time is then the average value 
'^■i^i ^ ''^{X'^'^j^). We deduce that the process defined by 

1=1 n=l \ jj^i I 
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is a local martingale. We finally get 



TV A\ 



i=l n=l 



N 



(6) 



which is exactly ([3]) (we recall that M and M implicitly depend on A^). 



Step 2: proof of inequality (|4| . 

The number of particles > 2 is still fixed. Let us first prove that, for any a S {1,2, . . .}, 
i'^^^n{'^tat^) ~ ^oi'^o) is a local martingale. In order to simplify the notations and 
since > 2 is fixed, we remove the superscripts N until the end of this step when there 
is no risk of confusion. 

The jump of the processes {At)t>o and A4 only occurs at the rebirth times, thus we 
deduce from (j6)) that 



/TV- 1 
N 



n=l 



N 



n=l 



Let us develop and compute each term — in the right hand side. For 

any n > 1, 



Tn Tn-) 



N - 1 



Tn 



+ Mr„-(V'r„-) 



N - 1 

N 



N - 1 

N 



On the one hand, denoting by i the index of the killed particle at time r„, we have 



Mrn-Mrn-, 



where 



N{N 



and, by the definition of M, 



1 



N -1 



N{N - 1) 



Tn- 



N{N -1)^^"-^^^--^ N 
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We then have 

On the other hand, we have 
A^- 1 



N-1 



+ 



N 



N-1 



N 

Finally we get 



jV - 1 \ ^^'^ 



1 



(Mr„-Mr^.)+Mr„-Mr„-, 



N- 1 



TV — 1 \ 



N - 1 



N 



N 



N- 1 



N 



We deduce that the process vtii^t) — J^o('0o) is the sum of two local martingales and we 
have 



N 



dMt + 



N 



dMt 



(7) 



Let us now bound both terms on the right-hand side of ([7| (where N is still fixed). In 
order to handle martingales instead of local martingales, we fix an integer a > 1 and we 
stop the particle system when the number of rebirths At reaches a, which is equivalent to 
stop the process at time (which is actually t^)- The processes M and M stopped at 
time Tq, are almost surely bounded by a||/||oo- By the optional stopping time theorem, 
we deduce that M and M stopped at time are martingales. On the one hand, the 
martingale jumps occuring at rebirth times A^r„ ~-^r„- are bounded by ||/||oo/-^ (since 
its jumps are the difference between two positive terms bounded by ||/||oo/-^) and the 
martingale is constant between rebirth times, thus 



E 



N - I fT/\Tc / jv - 1 ^ 



N 



N 



dMt 



< E 



< 



n=l 



N 



On the other hand, we have 




dMt 



N 



<E (MTArJ^ 



(8) 



1 ^ 



where 



m=0,n=0 
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If i 7^ j, then the expectation of the product of the martingales M*'" and M-'''" is 0, since 
the particles are independent between the rebirths and do not undergo simultaneous 
rebirths. Thus we have 

Assume now that i = j and fix m < n. By definition, we have 



which is measurable with respect to X^^^^^^^i then 

applying the optional sampling theorem with the martingale M^"^^ stopped at the uni- 
formly bounded stopping time T A Tq A t^. We deduce that 



\n=0 

<£('X;*TA,;W„,;)')+ll/ll?o 

\n=0 / 

<ll/llooi^ (E^ta.a(^ta.a)) + 



2 

oo • 



^n=0 



By ([5]), we have 



thus 



We finally have 



\n=0 / 
^((^TA.J')<2 




2 

oo ■ 



rfMt < (9) 



N ) I - N 

The decomposition ([7| and inequalities ([8| and ([9| lead us to 

U + V2)||/||, 



N 
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By assumption, the number of rebirths of the A^-particles system remains bounded up 
to time T almost surely. We deduce that T A converges to T when a goes to infinity 
(actually T At^ = T for a big enough almost surely). As a consequence, letting a go to 
infinity in the above inequality and using the dominated convergence theorem, we obtain 



Step 3: conclusion of the proof of Theorem [3 

By assumption, //^ {Pt ■) converges weakly to ^q{Pt.). Thus, for any bounded contin- 
uous function f : F ^ IR,_|_, the sequence of random variables (//^ {P^ Ip), fJ-Q {Pt f)) 
converges in distribution to the random variable (;Uo(-PrlF)> /Uo(-Pt/))- By (jlOp . we de- 
duce that the sequence of random variables {u^ {1f),i^t (/)) converges in distribution to 
the random variable {iiq{Pt'\-f)t l^o{PTf))- Now ^q{Pt'\-f) never vanishes almost surely, 
thus 



^f(lF) 

for any bounded continuous function / : F 
rem[T] 

We also deduce from (1101) that 



law 



MPrf) 



—7- 11+ , which implies the first part of Theo- 



E 



AT - 1\ 



N 



then, using (|10p and the triangle inequality. 



£;(|/x^(p^i^K(/)-M^(p^/) 



< 



1 + V2 



,2A ^ 2(l_+^/2)||/||, 



N 



Using the Cauchy-Schwartz inequality, we finally deduce that 



E 



f^^iP^f) 



which concludes the proof of Theorem [1] 



< 



E 



2(2 + ^/2)11/11, 



N 



□ 



3 Criterion ensuring the non-explosion of the number of 
rebirths 

In this section, we fix > 2 and consider A^-particles systems whose particles evolve 
as independent diffusion processes and are subject to rebirth after hitting a boundary 
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or after some exponential times. We begin by defining the diffusion processes which will 
drive the particles between the rebirths, then we define the jump measures giving the 
distribution of the rebirths positions of the particles. As explained in the introduction, 
the interacting particle system will be well defined if and only if there is no accumulation 
of rebirths in finite time almost surely. The main result of this section is a criterion 
(given by Hypotheses [1] and [2] below) which ensures that the number of rebirths remains 
finite in finite time almost surely, so that the A^-particles system is well defined at any 
time t > 0. 

For each i € {1, . . . ,N}, let Ei be an open subset of M'^* [di > 1) and Di be an open 
subset of M'^i (d^ > 1). Let {e\,Zl)t>o be a time-inhomogeneous environment- dependent 
diffusion process evolving in EiX Di, where G Ei denotes the state of the environment 
and Zj G Di denotes the actual position of the diffusion at time t. Each diffusion process 
{e\,Zf)t>Q will be used to define the dynamic of the i^^ particle of the system between its 
rebirths. By a time-inhomogeneous environment- dependent diffusion process, we mean 
that, for any i S {!,..., A^}, there exist four measurable functions 

Si : [0, + oo[y.Ei x Di ^ M'^' x M'^' 
rui : [0, + oo[xEi x Di ^ 
Ui : [0, + oo[xEi X Di^ M< x R'^'^ 
rji : [0, + co[xEi x Di ^ M'^'s 

such that (e*,Z*) is solution to the stochastic differential system 

del = Si{t,ei,Zi)df3l + m^{t4,Zi)dt, ej, e Ei, 
dZl = a,{t,ei,Zi)dBi + r,\t,e'„Zl)dt, Z^ G A, 

where the (/3*,i?*), i = 1, . . . ,N, are independent standard di -\-d[ Brownian motions. We 
assume that the process {e\,Zf)t>o is subject to two different kind of killings: it is killed 
when Zl hits dDi {hard killing) and with a rate of killing Ki{t,e\,Zl) > {soft killing), 
where 

Ki : [0, + oo[x£'i X Di H> M+ 

is a uniformly bounded measurable function. When a process (e*,Z*) is killed, it is 
immediately sent to a cemetery point d ^ EiX Di where it remains forever. In particular, 
the killed process is cadlag with values in Ei x DiU {d}. 

We assume that we're given two measurable functions 

5 : [0, + oo[xii^i X • • • X En x Di x ■ ■ ■ x Dn — ^ Mi{Ei x ■ ■ ■ x En x Di x ■ ■ ■ x Dn) 
and 

n : [0, + oo[xEi x---xEnx d{Di x • • • x Dn) Mi{Ei x ■ ■ ■ x En >^ Di x ■ ■ ■ x Dn), 
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where Mi{Ei x • • • x E^' x Di x • • • x Djv) denotes the space of probabihty measures on 
£^1 X • • • X Ejsf X Di X • • • X D^. These measures will be used to define the distribution 
of the rebirths locations of the particles, respectively for soft and hard killings. 

The A^-particles system is denoted by (Ot,Xt)j>o, where Ot = (o|, • • • ,o^) € EiX - ■ ■ xE]\[ 
is the vector of environments at time t and Xj = {X^, ■ ■ ■ ,X[^) ^ Di x ■ ■ ■ x Dn is the 
vector of positions of the particles at time t. In particular, the position of the i^^ particle 
of the system at time t > is given by (oJ,Xj) £ Ei x Di. 

Dynamic of the particle system 

The particles of the system evolve as independent copies of (e*,Z*), i = 1,...,N, until one 
of them is killed. Then 

- if two or more particles are simultaneously killed at this time, we stop the definition 
of the particle system and say that it undergoes a failure. This time is denoted by 

Tstop- 

- otherwise, 

- if the unique killed particle is softly killed, then the whole system jumps 
instantaneously with respect to the jump measure 5(t,Oj_,Xj_) 

- if the unique killed particle is hardly killed, then the whole system jumps 
instantaneously with respect to ^(t,Ot-,Xt_) 

and, in both cases, we say that the system undergoes a rebirth. 

After this operation and in the case of a rebirth, all the particles lie in F. Then they 
evolve as independent copies of (e*,Z*), i = 1,...,N, until one of them is killed (which 
leads to a failure or a rebirth) and so on. 

If the process does never fail, we set Tstop = +00. We denote the successive rebirth times 
of the particle system by 

n <T2 < ... <Tn< ... 

and we set Too = lim„_>.oo Tn. Since there is no natural way to define the particle system 
after time Tqo A Tstop, the A^-particles system is well defined at any time t > if and 
only if Too A Tgtop = +00 almost surely. In Theorem [2] below, we state that this happens 
if Hypotheses [1] and [2] below are fulfilled; this is the main result of this section. We 
emphasize that the above construction and the assumptions of Theorem |2] include the 
particular case of the Fleming- Viot type system studied in Section [2] 

The first assumption concerns the processes (e*,Z*), i = 1, . . . ,A^, which drive the particles 
between the rebirths. We denote by (pi the Euclidean distance to the boundary dDi, 
defined for all x G M"^* by 

(j)i{x) = inf \\x - y\\2, 
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where ||.||2 denotes the Euchdean norm on M^'i. For any a > and any i G {1, • • • ,N}, 
we define the boundary's neighbourhood by 

Df = {x £ Di, (piix) < a}. 

Hypothesis 1. We assume that there exist five positive constants gq, kg, cq and Cq such 
that 

1. (pi is of class on , with uniformly bounded derivatives, 

2. Ki is uniformly bounded on [0, + oo[x£'j x Di and Si,ai,mi and r]i are uniformly 
bounded on [0, + oo[xE'j x L)"", 

3. for any i G {1,...,A^}, there exist two measurable functions fi : [0, + oo[xi?j x Dl° — )■ 
M+ and Qi : [0, + oo[xEi x Df" M such that V(t,e,z) e [0, + oo[xEi x Df", 

X] = fi{t,e,z) +gi{t,e,z), (11) 

and such that 

(a) fi is of class in time and of class in environment/ space on [0,+oo[x£'j x 
-D"°, with uniformly bounded derivatives, 

(b) for all {t,e,z) G [0, + oo[x£;i x D^'' , 

\gi{t,e,z)\ < kg(j)i{z), 

(c) for all {t,e,z) e[0, + oo[xEiX , 

Co < fi{t,e,z) < Co and cq < fi{t,e,z) + gi{t,e,z) < Cq. 



Remark 5. We recall that the regularity of (pi near the boundary is equivalent to the 
(jk regularity of the boundary dDi itself, for all /c > 2 (see 0, Chapter 5, Section 4]). 

Remark 6. The euclidean distances (pi in Hypothesis [1] could be replaced by other func- 
tions with similar properties, in order to weaken the assumption. However, for sake of 



clarity, we do not develop this possibility and refer the reader to [l3| for an alternative 
class of distance like functions. 

Our second assumption is related to the rebirth measures 7i and S. We emphasize that 
if a rebirth due to hard killing happen, then at most one particle hits its associated 
boundary dDi. This implies that, at a rebirth time, the whole set of particles hits one 
and only one of the sets Di, i = 1, ■ ■ ■ ,N , defined by 

A = {(^^i, • • • ,XAr) e d{Di X ... X Dn), Xi G dDi and Xj G -Dj,Vj / i] . 

With this definition, it is clear that the i*^ particle undergoes a hard killing leading to a 
rebirth if and only if (•,0,X) hits [0, + oo[x (i?i x . . . x E^) x "Dj. In particular, the values 
taken by % outside these sets does not influence the behaviour of the particle system. 
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Hypothesis 2. We assume that 

1. There exists a non- decreasing continuous function h : — )■ M+ vanishing only at 
such that, Vi G {1,...,N}, 

inf ^(t,e,2;i,...,2;Ar) x ... x E^at) x Ai) > po, 

(t,e,(xi,...,xjv))e[0,+oo[x(_EiX...x_Bjv)xI', 
where pQ > is a positive constant and Ai C Di x . . . x Djv i/ie set defined by 
Ai = {(yi, . . . ,yN) I 3i / i such that (j)i{yi) > h{<j)j{yj))} . 

2. We have 

inf n{t,e,xi,...,XN)iBe,x) = 1, 

(i ,e,x) 6 [0,+oo [ X (El X ■ • • X _Bjv ) X 



where B^^x C (E'l x . . . x En) x (Z^i x . . . x L)jv) is defined by 
Be,x = I ie,x')\\/j, (j)jix'j) > (pjixj) 



\l fj{t-,GjiXj) 



fj being the function of Hypothesis [Jl extended outside [0, +oo[xEj x Dj° by the 
value cq. 

Let us explain Hypothesis [2] 

- The set Ai is a subset of Di x ... x Djv such that if the i*'* component of an element 
(?/!)■■■ )2/n) G is near the boundary dDi (that is if (j)i{yi) <^ 1), then at least one 
another component yj fulfils h{(f)j{yj)) <C 1 and thus (pDiVj) ^ 1- This implies that if 
the rebirth position of the i*'* particle after a hard killing is located near the boundary 
dDi, then, with a probability lowered by po, at least one other particle is located near 
the boundary. 

- The set i?e,x is a subset oi EiX . . .x E^ x Di x . . . x such that for any {e' ,x') G Be^x, 
the component x'j lies farther from the boundary than Xj, for any j G {1, • ■ • This 
means that after a hard killing, the rebirth position of any particle lies farther from the 
boundary after the rebirth than before the rebirth. 



Remark 7. Hypothesis [2] includes the case studied by Grigorescu and Kang in [13| (with- 
out environment dependency or time inhomogeneity) , where 

n = '^Pijixi)5xj, V(xi,...,XAr) G A- 

with '^j^iPij{xi) = 1 and infjgji /^r} .;,.gg£)Pjj(xi) > 0. Indeed, in that case, the 
rebirth position of any killed particle is the position of another particle. It implies that 
Hypothesis |2] is fulfilled with pQ = 1 and h{(f)j) = (pj. This case also includes the Fleming- 
Viot type system of the approximation method proved in Section [2j 

We're now able to state the main result of this section. 
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Theorem 2. If Hypotheses[I\and\^are satisfied, then Tstop^Too = +00 almost surely. In 
particular, the N -particles system with rebirth (•,0,X) is well defined at any time t >0. 

Remark 8. An other interesting example of diffusion processes with kilUng is given by 



reflected ones (see for instance |14| and |25|). Thus a natural question is whether our 
non-explosion result holds true when the diffusion process (e*,Z') is reflected on dDi 
and killed when its local time on dDi reaches the value of an independent exponential 
random variable. In fact, the only difference lies in the Ito calculus developed in the end 
of this section, where a local time term should appear due to reflection. The generality 
of Proposition [3] (which concludes the proof) allows to handle this additional term. As 
a consequence, our proof would remain perfectly valid under the setting of reflected 
diffusion processes. 

Proof of Theorem Since the rate of killing Ki is assumed to be uniformly bounded for 
all i € {!,..., A^}, there is no accumulation of soft killing events almost surely. As a 
consequence, we only have to show the non-accumulation of hard killing events and can 
assume up to the end of the proof that = for all i E {!,..., A^}. 

In Subsection 13. 1| we prove that, conditionally to Tgtop A Tqo < +00, there exists i 7^ 
j G {1, . . . ,A^} such that the pair of semi-martingales {<j)i{Xl),<j)j{Xf)) converges to (0,0) 
when t goes to Tgtop A Too- This leads us to 



P {Tstop A Too < +00) < 5Z ^ (^O'' < 

l<i<j<N 

where Tq is the stopping time defined by 

= inf \t G [0, + oo[, 3sn > t such that lim {^KJ^^jiXij) = (0,0)| , 

where (s„) runs over the set of non-decreasing sequences. 

In a second step, we conclude the proof of Theorem |2] by proving that, for any i 7^ j, 

P{T^^ < +00) = 0. 

The proof of this assertion is itself divided into two parts: in Subsection 13. 2| we com- 
pute in details the Ito's decomposition of <j)i{X^) and (j)j{X^); in Section HI we prove 
that the just obtained Ito's decomposition implies the non-attainability of (0,0) for 
{(f)i{X^),(j)j[X^))^ that is Tq = -\-oo almost surely, which concludes the proof. 

Remark 9. The non-attainability result used in the proof has its own interest and is thus 
stated as Proposition [3] in the independent Section H] 
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3.1 Explosion of the number of rebirths implies simultaneous conver- 
gence to the boundary 

Let us prove that, under Hypotheses [T] and |2] 

P {Tstop A Too < +00) < ^ {^0 < +00) . (12) 

l<«<j<7V 



If Tstop < +00, then at least two particles X* and hit their respective boundaries at 
time Tstop (by definition of this stopping time). It impUes that (/>j(X*^^^^_) = (pji^r^^^p-) = 
0. We deduce that 

{Tstop < +00} C U {^0^ < +°°} • (13) 

l<i<j<Af 

Define the event £ = {tqo < +00 and Tstop = +00}. It remains us to prove that, up to a 
negligible event, 

£c u {n'<+^}- (14) 

l<i<j<N 

Conditionally to £, the total number of rebirths converges to +00 up to time Too- Since 
there is only a finite number of particles, at least one of them, say io, is killed infinitely 
many times up to time Tqo ■ We denote by {T^)n the sequence of rebirth times of the 
particle zq (we thus have r^" < Too < 00 for all n). In a first step, we prove that (f)i^^{X^\^) 

converges to when n — )• 00. In a second step, we prove that there exists a random index 
io 7^ ^0 such that (j)j^^{X-'°^) converges to when n — >■ 00. This will imply (|14p and thus 

Inequality ([T2|). 

Step 1: convergence of (j)ig{X^1^) to when n — )• 00. 
We prove that, up to a negligible event, 

£cU,{XX)^0}. (15) 

Ito's formula will be a very useful tool in our proof. However, since (pi^ is of class 
only on the set -D^^", one can not apply directly the Ito's formula to </)j(, (X*" ) . In order 
to overcome this difficulty, let 7 : ]R_|. 1— t- IR+ be a function such that 

7 is non-decreasing with uniformly bounded derivatives, 

-f{u) = ^, Vn G [ao, + oo[. 

The function 7 o is of class on the whole set Di^ with derivatives 

^W^(^) = I 7'{Mx)f-^{x), if X G Dl<; 
dxi I otherwise. 
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and 



dxidxk 



7"(</'.o(^))S^(-) +y(<A.o(x))£fe(x), if X G Dl 



I 0, otherwise, 

where k,l £ {1, . . . ,d^^}. By definition of 7 and by Hypothesis [1] these derivatives are 
uniformly bounded over Di^ and one can apply the Ito's formula to 7 o 

Applying Ito's formula to 7 o <j)i^(X!°) between times and t^'^d we obtain 



1=1 



dxi 



dxidxi 



+ 



1 " 

i,k=i -'^^ 
1=1 ''^^ 

i[K{xT))-i[K{m)- 



dx, 



dt 



n + 1 



By Hypothesis [21 (j)i^^{X'^°) can only have non-negative jumps. Since 7 is non-decreasing, 
we deduce that the last sum is non-negative. Since r^'^i — converges to and the 
derivatives of 7 o (f)^^ are uniformly bounded, we also deduce that the integrals in the 
right term converge to when n goes to -|-oo. We thus have, conditionally to the event 



liminf7U,„(X\ 



^ ' n 



> 0. 



(16) 



But T^^i is by definition a rebirth time for X*o and thus {X^l^ ) = (we recall that 

we assumed without loss of generality that Kj = 0, so that all rebirths are due to hard 
killing). Since the continuous function 7 is equal to at 0, we deduce that 



7 { c/>.o (^^^0 ) ) - 7 ( '^.0 ))=-l[<P^o iX% ) ) < 0. 

'''n + l ^ 



This and the limit (|16p imply that 

lim J (4>^,{X%[ 
which implies (|15p . 

Step 2: there exists jo 7^ iq such that (j)j^^{X-"-l^) converges to 0. 
We prove now that, up to a negligible event, 

e C |3jo / io, K{X\) > h{(t>j,{xX)), for infinitely many n > l| , 



(17) 
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where h is the continuous function of Hypothesis [2] We denote by An the event 

A = < +00 and, Vj / io, 4'^o{X\) < h{^j{X\J)} . 
We have, for all 1 < A; < /, 

By construction of the rebirth mechanism of the interacting particle system and by the 
first point of Hypothesis [2l 

< 1 - Po, 

where and ^-re defined in Hypothesis |2] By induction on I, we get 



p(r\An] <(1 



Po)^"^ VI < A: < L 



Since po > 0, it yields that 

p(u n-^n) =0- 

\k>ln=k J 

As a consequence, if £ happens, then, for infinitely many rebirth times r^", one can find 
a particle j 7^ such that > h{(j)j{X-' ^^J). Since there is only a finite number 

of other particles, one can find a particle, say jo (which is a random index), such that 

£ C Ui,{X\) > h{^j,{X\)), for infinitely many n > l| . 
Let us now conclude the proof of (|14p . We deduce from (|15p and Step 2 that 

<5 C ^(<?^7n (^"'in )) ^ 0; to a subsequence > , 

C \ (/)j,AX-^° ) 7- 0, up to a subsequence > . 

Using again (|15p . we obtain that, up to a negligible event, 

£ C \3jo / io, lim ((^i,j(X*'^J,(?:>j3(X^°J) = (0,0) up to a subsequence! 

N 

c u {n' < +00} . 

This and (|13p clearly implies (|12p . 
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3.2 Non-convergence to (0,0) for {(j)i{X'),(f)j{X')) 



Fix i ^ j & {1,...,A^}. In this subsection, we prove that 



P ro*^<+oo =0, (18) 



which will conclude the proof of Theorem |2] 

Let us introduce the sequence of stopping times (tn)neN ) defined by 



to = inf{t G [0,T,top A Too[, ^J U^iy + <PjiXiV < ao/2} 
and, for all n > 0, 

t2n+l = inf{t G [t2n,Tstop A Too [, \J (pi{Xl)'^ + (/)j(X/)2 > Oq} 



t2n+2 = mf{t G [t2n+l,Tstop A Too[, ^ (Ai(^l)^ + <Pj{XlY < ao/2}. 

By definition of the sequence (tn)n>i) we clearly have 

f 0i( X^') < ap and (/>j(^ t ) < ao, Vt G Un.>o[*2n,i2n+i[) 
)i(X*)2 + (l)j{xiY > ao/2 otherwise. 

In particular, {(j)i{Xl) (XD) cannot converge to (0,0) during the interval of time t G 
[t2n+i ,t2n+2[ and thus Tq ^ [t2n+i,t2n+2[ almost surely, for any n > 0. We deduce that 

+ 00 

P {t^' < +CXd) <Y,P (t^' G [t2n,t2n+l[) + P (t^'' > ^n) • 

n=0 

Our first step consists in proving that the sequence (tn)n>o converges to +oo almost 
surely. In a second step, we prove that P ( Tq-' G [t2n>i2n+i[) = for each n > 0. 



Step 1: {tn)n>o converges to +oo almost surely. 
Since t„ is non-decreasing, it is clear that 

{{tn)n>o does not converge to + oo} C < supt„ < +00 

[n>0 

C \ t2n+l — t2n > > . 

Using the same Ito's decomposition as in the first step of Subsection 13. 1| we obtain that, 
up to a negligible set, 

jWi - t2n > 0| C jliminf 7 o (j^kiX^ ) -70 <j)k{XU < o| for k = 
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Since 7 is non-decreasing, we deduce that, up to a negligible set, 

IWi - t2n > 0} C (liminf ) - MXtJ < 0| for k = 

We thus have 



c |iimmf ^jMxL„r + ■PAK_,.„r - \/MxLr + -PAK.)' £ »| ■ 

But the right continuity of the particle system, the continuity of (pi and and the 
definition of the sequence imply that 



supt„ < +00 

n>l 



Finally, we deduce that {(tn)n>o does not converge to + 00} is included in a negligible 
set, so that {tn)n>i converges to +00 almost surely. 



Step 2: we have P [Tq G \t2n:'t2n+iVj — f^'^ '^^y n > 0. 
Fix n > 0. We define the positive semi-martingale y* by 



OO V 1^ if t > t2n+l - t2n- 



(19) 



We also define the local martingale and the adapted process ¥ by 



» ft/\{t2n + l-t2n) 



1=1 



dxi 



t2n+sl I ' 



and 



- ltG[0,t2„ + l-t2n[ 



Z,fc=l ' 



With these definitions, we deduce from the Ito's formula that 

dYl = dMi + h\dt + y/ - y/\ 



(20) 
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We define the adapted processes vr* and /?* by 

I fiihn + t,ol^^_^_^,Xl^^_^_^), if t < t2n+i - hn, 

I Co, if i > t2n+l — t2n 



Pi 



and 

5i(*2n + t,0^2^+j,Xj*2„_^j), if t < t2n+l " *2n, 
^— Co, if i > t2n+l — ^271, 

where fi and are given by Hypothesis [1] An immediate computation leads us to 

d{A'P)t = (4 + pi)dt, Vt > 0. (21) 

The process vr* is a semi-martingale. We deduce from the Ito's formula and from the 
regularity of fi over [0,-|-oo[x£^j x L)^" that there exist a local martingale A^* and a finite 
variational process U such that, for all t >0, 



diri = dNl + dL\ + t^I - 7rj_. 
We define the adapted process ^* by 



(22) 



di 



Lte[0,t2n + l-i2n[ 



E 



9et. dei 



[SiS*]kl + X] 



dxk dxi 



where each function involved is implicitly evaluated at (i2n + t,ol^^_^f,Xl^^_^,f) . We thus 
obtain 

{N')t = Ctdt. (23) 
Similarly, we define the processes , , V , , p> , , U and 
By Hypothesis [2] and by the definition of y* and , it is clear that 



y: 



— ^ > and * 



>0. 



(24) 



Moreover, the uniform bound assumptions and the regularity assumptions of Hypothesis [T] 
imply that the processes ¥, V, vr*, ir^ , p^, , ^* and are uniformly bounded. In 
particular, there exist some positive constants 6oo, and such that, for any t >0, 



(25) 



bl > -boo and > -b^o, 
vrj < ttI + Ip^I < and 4 < vr^' + \pl\ < C^, 
[ < Q and C{ < C^. 

Setting = cq and ko = kg \/ co/ao, we deduce from the fourth point of Hypothesis [T] 
and the definition of vr', vr-', and p^ that, for any t > 0, 



Ctt < Ti'f and < tt^, 
|pj|<A;oy/and|p^|<A:oy/. 



(26) 
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By the independence of the particles between the rebirths, we also deduce that, for any 
t > 0, 

{M\M^)t = 0. (27) 



We claim now that the decompositions 1^0^, ([21]), (l22|) and (l23|) of , (M'), (M^), 
7r\ vr^ {N') and {N^), together with the properties (l24|), ([25]), (|26]) and 1^, imply that 
(y^,y^) never converges to (0,0) almost surely. This is proved in the next section, where 
a general criterion for non-attainability of (0,0) for non-negative semi-martingales with 
positive jumps is stated (see Proposition |3]). 

Finally, we deduce that Tq ^ [t2n,t2n+i[ almost surely, for all n > 0, concluding the 
proof of Theorem [2] 

□ 



4 Non- attainability of (0,0) for semi- martingales 

Fix T > and let (Xt)te[o,T]j ^ — 1)2, be two uniformly bounded non-negative one- 
dimensional semi-martingales. This means that there exists a constant yoo > such 
that 

< y/ < Voo, Vt > and z G {1,2} almost surely (28) 

and that there exist a continuous local martingale (M|)jg[o,T] ^^'^ ^ continuous finite 
variational process I such that 

dV; = dMl + dli + y/ - y/', Vt > O and i E {1,2} almost surely. (29) 

In this section, we give a sufficient condition for (y^,y^) not to converge to (0,0) up to 
time T. More precisely, we set 

To = inf [t G [0, + oo[, 3s„ > t such that lim {Y} ) = (0,0)1 , 

where (s„) runs over the set of non-decreasing sequences, and give a sufficient criterion 
for To = -|-oo almost surely. 

Criteria for non-attainability of (0,0) already exist for time homogeneous stochastic dif- 
ferential equations (we refer the reader to the papers of Friedman flo| . Ramasubrama- 



nian 



2l| and Delarue 0). esult is a generalization, since we do not assume at all 

that y* is given by a stochastic differential equation, nor that it has any time-homogeneity 
property. In particular, it is not required for y* to be a Markov process, which is of first 
importance in view of our main application (the processes (j)i{Xl) in Section 3 do not 
fulfil the Markov property). Our proofs are inspired by the recent work of Delarue Q, 
who obtains lower and higher bound for the hitting time of a corner for a diffusion driven 
by a time homogeneous stochastic differential equations reflected in the square. 

Here is our main assumption. 
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Hypothesis 3. For each i = 1,2, there exist an adapted process 6* and a non- decreasing 
adapted process Kl such that 

dlf = bidt + dKf. 

Moreover, there exists a non-negative semi-martingale vr*, whose decomposition is 

d7:l=dNi + dLl + 7rl-7rl, 

where is a continuous local martingale and U is a continuous finite variational adapted 
process, and there exist two adapted processes p\ and Q and some positive constants 
6ooifcoiC7r,C7r, such that, almost surely, 

1. d (M^)^ = {-kI + p\)dt and d {N')^ = Qdt, 

2. c^<Til + p\< C^, \p\\ < koY^', < Q and b\ > -boo for all t G [0,T] 

3. (^M^,M'^'j is a non-increasing process. 
4- we assume that, for any time t > 0, 



> 0. 



Remark 10. The third point of Hypothesis [3] has the following geometrical interpretation: 
when an increment of is non-positive (that is when goes closer to 0), the increment 
of is non-negative (so that goes farther from 0), as a consequence (M"'^,M^) 
remains away from 0. A nice graphic representation of a very similar phenomenon is 
given by Delarue's 0, Figure 1]. 

We define the adapted process ^ by 

(30) 



2 =v < , 

For all e > 0, we also define the stopping time = inf{t G [0,T], $t > e~^}. Since ttj is 
uniformly bounded below by Ctt, goes to infinity when {Y^ ,Yi) goes to (0,0), so that 

To = lim almost surely. 

e-5>0 

We are now able to state our non-attainability result. 

Proposition 3. If Hypothesis is fulfilled, then Tq = -\-oo almost surely. In particular, 
(y^,y^) doesn't converge to (0,0) in finite time almost surely. Moreover, there exists a 
positive constant C which only depends 071 ^ko^Cj^ jC'j^ ^ such thcitj for all £ ^ ^ 
and any stopping time T, 

P (T, < T) < pr^^ {Ei\L\ + + T)) , 
where \L^\t is the total variation of U at time T and $o is defined in (|30p . 
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Proof of Proposition Since converges to Tq when e — )• 0, if the last part of Proposi- 
tion |3] is fulfilled, then P(To < T) = for any deterministic time T and thus Tq = +00 
almost surely. As a consequence, we only have to prove the second part of the proposition, 
for a given stopping time T and a fixed value of e > 0. 

The proof is divided into several steps and is organised as follows. We assume that 
(M^,M^) = in the three first steps and consider the general case in the last one. In 
Step 1, we compute in detail the Ito's decomposition of the semi-martingale In Step 2, 
we introduce the function F : R 1— t- E, defined by 



/•r 

F{r) = / exp {Cpe-') ds, 
Jo 



where Ci? > is a positive constant that will be fixed later in the proof. Using Ito's 
formula, we prove that there exists a local martingale H and a positive constant C > 
such that 

F{^t) - i^(^o) <Ht + C {\L^\t + \L^\t + t) . (31) 

In Step 3, we conclude the proof of Proposition [3] in the particular case {M^,M'^) = 0. 
In Step 4, we prove that the result remains true in the general case. 

Step 1. Ito's decomposition of the semi-martingale {^t)t&[o^Te[- 
Let us introduce the function F defined by 

r : r:^ X e:; X R+ X E+ ^ r 

1 f 

{ai,a2,xi,X2) i-> -- log — — 

2 \ai 02 

so that = r(7r^,7r|,y^-'^,y^^). The successive derivatives of the function F are 

dY 1 _2 2 2r d'^Y _3 2 2r ^ —4 4 4r 
^=2"* ,^ = -a, x,e +-a, x^e , 

^ —2 2r —3 3 4r t/ J- _2 2 4r 

Q^,=ai XiC -a, x^e > = -«i «j , 

9 ^ ~ -1„-1 4r ^ r _ ^ „-2„-2„2_2 4r 



In particular, we have 

rP-V 

Vc -L , , > 1 op 2 2 4r 

— ^(ai,a2,a^i,a^2)ai = 2^ e 2q • x^e 



i=l,2 i=l,2 



a\ a2 
.2e2r + 26^^6-2^ = 0, 
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so that, for any t G [0,Te[, 

XI -^('^t,'^tXt^Xt^)'^l = 0) almost surely. 



^=l,2 



Using the previous equalities and the Ito's formula, we get 



i=l,2 ^* i=l,2 i=l,2 ^* 



and 



i=l,2 i=l,2 ^^^i> 

+ -y f--e2** + 2plje4**l pjdt + - y ^^e^*M(M\M^)i 
j=i,2 ^ * ^ *^ ^ ■i^je{l,2} "i"t 



(32) 



i=l,2 i=l,2 



i=i^2'^^"t> i^ie{i,2} ^^tv^t'' 



^ie;? 5; proof of . 

Let Cp > be a positive constant that will be fixed later in the proof and define the 
function F : E i-)- E by 

F{r) = [ exp {Cpe-') ds. 
Jo 

Setting Dp = Ci?yoo\/2/\/c^) one can easily check that 

1 < F'(r) < e^^ and F"(r) = -Cpe-'' F' {r), Vr > --log (l^^ , 

2 V Cjr I 
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where — ^log ^2^^ is a lower bound for the process By Ito's formula, we deduce 

that, for any t S [0,re[, 

F{^t)-F{^o)= f F'{^s)d^';-^ /* e-*=F'($,)d V F{^,)-F{^,.), (33) 

•^0 ^ -^0 0<s<t 

where is the continuous part of d^s- Our aim is to prove that 

F{^t) - F(^o) <Ht + C {\L\ + \L\ + t) , 

where H is the local martingale defined for any t > by 

= - E / ^e''''F'{^s)dM: + ^ / lh^e''^^F'{^,)dNl. 



1=1,2 



This is done below using (|32p and proving lower or higher bounds for each term on the 
right hand side of (|33p . respectively in step 2a, in step 2b and in step 2c. 

Step 2a. Let us prove that there exists a positive constant C' > 0, which does not depend 
on Ci?, such that 

/ F'{<^,)d^l <Ht + C' e''^F'{^,)ds + ^ {\L'\t + \L\) + ^^-P^. (34) 

Jo Jo ^Ctt ^C^ 

Since is non-decreasing, we have 

- E ^^e"''F'{^s)dKl<0. 

One can easily check that l^*e** < \/vr^, then ^e*' < Since h\ > — 600 we have 



Since the derivative F takes its values in [l,e^l and since , \L e ' < — , we deduce 



that 



By the second point of Hypothesis [3l we have \p\\e^^ < koY^e'^' < /coy^ < /co^/CV- As 
a consequence, using that (^^p ^ ^ ~ ^rj ^ C71-, we deduce that 
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By the Kunita-Watanabe inequality (see [22|, Corollary 1.16 of Chapter IV]) and by 
Hypothesis |3] (first and second points), we get, for all predictable process hs and any 
couple i,j G {1,2}, 



and, similarly. 



< 



^ j\hs\{M^),^ j\hs\{W), < yCVQ \hs\ds, 



hs{M\M^)^ 



< / \hs\ds and 



< Q / \hs\ds, 



In particular, for any t G [0,re[, using Vfyj-e " < — and vr* > c^, we deduce that 



1=1,2 



2(vr: 



i\4 



and 



2(vr|)2(vr^2 

We also deduce that 



j=l,2 



s g2$. 



j\3 



F'(a>,)d(A'r,iV^)^ < 



2yc;Q /•* 



3/2 



Jo 



and 



4 E 



Setting C" 



vr*(7r^) 
26oo I 3fcoVCV , 3^C^C; 



e^*=F'($,)d(Ar,iV-'\ < 



3/2 



+ 



+ 



3/2 



> and using equation (|32p (recall that we 



assumed (M^,M^) = 0), we deduce that the higher bound (|34p holds almost surely, for 
any t G [0,r,[. 

Step 2b. We prove now the following lower bound for e^^" F' {^s)d {^) g, 

1^ e--^F'i^g)d {^)g > ^ 1^ e-^F'i^g)ds- (^^ + ^ + '-^^^ e^n. (35) 
One can easily check (using the second point of Hypothesis [3| that 



(Yi\2 „2$s , 
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Using that vr* > Ctt, we deduce that 

ft (vi\2 



The process (A^*) being non-decreasing and F' being positive, we have 



^e**^F'($,)e-*=d(Af*)^ >0. 
The same apphcation of the Kunita-Watanabe inequahty as above leads us to 



2c2 



since the definition of ^ imphes that it is uniformly bounded below by 0. We also have 
and 



- y f iML^'^^spi^^ y-^s^/M\N^) >-^[^^e^n. 

We finally deduce by ([32]) (where (A/\Af2) = 0) that ^ holds. 
5^6^ 2c. The jumps of are non-positive and F is increasing, thus 

^ - F($,.) < 0. (36) 

0<s<t 

Finally, by dMI), (I35I) and dM]), we deduce from <^ that 
m) - F($o) < + (C' - ^) e*»F'($,)d. + ^ (IL^I, + 



+ 



' (3 + c^)Cg fepc^yc; c^7cvq \ ^ 

24 2c. + 3/2 1^ 



Choosing Cp = 2Ct^C' /c-,^^ we have proved that there exists C > such that, for any 

t e [o,r,[, 

F($f) - F($o) <Ht + C {\L\ + \L\ + t) . 
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Step 3: conclusion in the case (|M^,M^^ = 0. 

Let {9n)neK be an increasing sequence of stopping times which converge to +00 such that 
(i7t)fg[o,0^] is a martingale for any n > 0. Let {O'^)n>o be the non-decreasing sequence of 
stopping times defined by 0^ = inf{i G [0,r], 3i st d\U\s > n} AT; since U is a finite 
variational process, {6'^) converges to +00 when n goes to 00. Using Step 2, we deduce 
that ^ 

^ ( / exp(e^^e^")du ) < C7E(|Li|e„ +\L\, +T). (37) 

Remark that ^tAT Ae' ab" reaches if and only if < 6' , A 6'',,. Using the right 
continuity of , Y'^, vr^ and vr^, we thus deduce that 

P(T,<a;,AC') = P{^T.Ae'^,Ae';„-^o>e-'-^o) 

( /■*T,Ae',Ae"„ „ \ 

< P\J " " exp(e^^e-")dn > e"^ - ^-oj , 

since r — q < JJ' exp(e*^^e~")(iu for all < 5 < r. Finally, using the Markov inequality 
and (|37p . we get, for all > ^0, 

P (T, < 6'^, A < -T^^Ce{\L\,^ + \L\,, + r) . 

Since (0^/) and (^"//) converge to +00 almost surely, letting n' and n" go to 00 implies 
the last part of Proposition |3l which implies the whole proposition. 

Step 4-' conclusion of the proof in the general case. 

Assume now that {^iVf^,M^^ is non-increasing. We define $J as the process starting 
from $0 and whose increments are defined by the right term of (|32p but removing the 
d{]VP,]VP) terms. We also define T^' = mi{t > 0, $j > e~^}. On the one hand, the same 
calculation as above leads to 

P{T^<T)< -—^-—CE(\L^\t + \L^\t + T). (38) 
e ^ - $0 

On the other hand. 

Since (M^,M^) is assumed to be non-increasing (third point of Hypothesis [S]) , we deduce 
that ^ ^'f It yields that T^' < almost surely, so that Proposition |3] holds even if 

□ 
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